Magnetism, superconductivity, and pairing symmetry in Fe-based superconductors 
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We analyze antiferromagnetism and superconductivity in novel Fe— based superconductors within 
the itinerant model of small electron and hole pockets near (0,0) and (tTjTt). We argue that the 
effective interactions in both channels logarithmically flow towards the same values at low energies, 
i.e., antiferromagnetism and superconductivity must be treated on equal footings. The magnetic 
instability comes first for equal sizes of the two pockets, but looses to superconductivity upon doping. 
The superconducting gap has no nodes, but changes sign between the two Fermi surfaces (extended 
s-wave symmetry). We argue that the T dependencies of the spin susceptibility and NMR relaxation 
rate for such state are exponential only at very low T, and can be well fitted by power-laws over a 
wide T range below Tc. 

PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.25. Ha 



INTRODUCTION 

Recent discovery of superconductivity in the iron- 
based layered pnictides with Tc ranging between 26 and 
52K generated enormous interest in the physics of these 
materials [l|, l^, y, U, l5[- The superconductivity has 
been discovered in oxygen containing RFcAsO (R=La, 
Nd, Sm) as well as in oxygen free AFe2As2 (A=Ba, Sr, 
Ca). Like the cuprates, the pnictides are highly two- 
dimensional, their parent material shows antifcrromag- 
netic long-range order below 150K [l|, [y, 0, H, Q, and 
superconductivity occurs upon doping of either electrons 
[l, 2, 3, 4] or holes [5| into the FcAs layers. 

The close proximity of antiferromagnetism and super- 
conductivity fueled early speculations that the physics 
of the pnictides is similar to the cuprates, and involves 
insulating behavior [lO|, [llj, ll2|. However, there is a 
growing consensus among researchers that Mott physics 
does not play a significant role for the iron pnictides, 
which remain itinerant for all doping levels, including 
parent compounds, in which magnetic order is of spin- 
density-wave (SDW) type rather than Heisenberg anti- 




ferromagnetism of localized spins [13|, [l^. This is ev- 
idenced by, e.g., a relatively small value of the ob- 
served magnetic moment per Fe atom, which is around 
12 — 16% of 2yLtB 0, Isl- In another distinction to the 
cuprates, electronic structure proposed by band struc- 
ture calculations [l5|, ll6|, |l7|, [l^, ll9| and supported by 
ARPES [20|, |21| consists of two small hole pockets cen- 
tered around F point (p = (0,0)) and two small electron 
pockets centered around M point (p = Q = (tt, tt)) in 
the folded Brillouin zone (BZ) (two Fe artoms in the 
unit cell, we set interatomic spacing a = 1) 

In this paper, we address three issues for the pnic- 
tides: (i) what interactions cause SDW order and su- 
perconductivity, (ii) what is the gap symmetry, and (iii) 



FIG. 1: (color online) A simplified FS geometry of doped Fe- 
based superconductors, used in the present work. At zero dop- 
ing, the Fermi surface consists of an electron pocket around 
(7r,7r) (black solid curve), and a hole pocket of roughly equal 
size around (0,0) (blue solid curve). We neglect in this work 
the fact that there are two hole and two electron pockets. In 
this system, there is a near-perfect nesting between hole and 
electron pockets (moving a hole FS by (tt, tt) one obtains a 
near-perfect match with an electron FS). Upon electron dop- 
ing, the size of the electron pocket increases (dashed blue -^ 
black), and this breaks the nesting. -|-A and —A are the 
values of the superconducting gaps on the two FS for s^ su- 
perconducting state. 



what are the implications of the gap symmetry for the ex- 
periments in the superconducting (SC) state. We argue 
that both magnetic and pairing instabilities are deter- 
mined by the same interband pair hopping which trans- 
forms two fermions near the hole Fermi surface into two 
fcrmions near the electron Fermi surface (and vice versa). 
This interaction may be weak at high energies (of order 
bandwidth), but it flows under rcnorm-group (RG) and 



ultimately determines the couplings in both SDW and 
Cooper channels at low-energies. When electron and hole 
pockets are nearly identical, SDW instability occurs at a 
higher T. When the near-identity is broken by either 
hole or electron dopings, the Cooper instability comes 
first. This pairing interaction sets the gaps in hole and 
electron pockets to be of equal magnitude A, but of op- 
posite signs (an extended s— wave symmetry, s+). The 
ratio 2A/rc = 3.53 is, however, the same as in BCS the- 
ory, as there is no angular variations of the gap along the 
FS. 

A fingerprint of s^ gap symmetry and near-equal elec- 
tron and hole pockets is the existence of a magnetic col- 
lective mode inside the gap for momenta near Q (a spin 
resonance), whose dispersion (^o + (w|^/2)(q— Q)^)"^^^, 
where vp is the Fermi velocity, closely resembles 
Anderson-Bogolubov mode in uncharged superconduc- 
tors. Another fingerprint is a strong reduction of the 
nuclear magnetic resonance (NMR) relaxation rate l/Ti 
in the clean limit, due to vanishing of the coherence factor 
for X {^.^^)/^ for q = Q. We argue that in this situa- 
tion, 1/Ti is predominantly due to impurities, which are 
partly pair-breaking even when non-magnetic. We show 
that, in the presence of impurity scattering, l/Ti is ex- 
ponential in T only for very low temperatures, and over 
a wide range of T < Tc is well described by 1/Ti ex T^, 
as if the gap had nodes. Over the same range of T, the 
uniform susceptibility is near-linear in T. 

Our results partly agree and partly disagree with some 
earlier works on Fe-pnictides. Mazin et al. 18| and 
Gorkov and Barzykin [22| conjectured that the pairing 
symmetry should be s+. Our results ag ree with theirs 
and also with Eremin and Korshunov [23[ , who analyzed 
numerically the magnetic response at Q within RPA for 
an s+ superconductor and found the resonance peak be- 
low 2A. Cvetkovic and Tesanovic [2j| noticed that for 
identical electron and hole pockets. Cooper and particle- 
hole channels become undistinguishable and should be 
treated equally - the notion we share. Wang et al [25|,|26| 
performed numerical RG study of the pairing symmetry 
and found an s+ gap symmetry for two-band model and a 
conventional s— wave symmetry for five-band model. We 



can only compare the results for the two-band model, for 
which we also found an attraction in in s"*" channel. There 
is, however, an important difference between our results 
and those of Wang et al. In our case, the bare interac- 
tion in s+ channel is repulsive, and attraction emerges 
only below some energy scale, due to RG flow of the 
coupling. In their analysis, the bare interaction is zero, 
and attraction emerges already after an infinitesimal RG 
transformation. Lorenziana et al. [27| used unrestricted 
Hartree-Fock approximation and studied possible phases 
that may compete with superconductivity in FcAs layers. 
We found that SDW is the main competitor, but CDW 
with complex order parameter is close second. 

On experimental side, ARPES(28l.[29l.[30t and Andreev 
spectroscopy[31j measurements have been interpreted as 
evidence for a nodeless gap, while NMR data were argued 
to follow 1/Ti ex T^ and were interpreted as evidence for 
a d— wave gap [32, l33[ or multiple gaps 3J|. Our results 
show that the T dependence of 1/Ti in a dirty s"*" super- 
conductor mimics T^ over a wide range of T and become 
exponential only at very low temperatures. 



THE MODEL 

We model iron pnictides by an itinerant electron sys- 
tem with two electronic orbitals, and we assume that the 
hybridization between the orbitals leads to small hole 
and electron pockets located near (0,0) and (7r,7r), re- 
spectively in the folded BZ (two Fe atoms per unit cell) 
(FiglT)). The extension to a more realistic case of four (or 
even five) orbitals and two hole and two electron pockets 
is straightforward, and docs not lead to new physics ex- 
cept for a magnetically ordered state, where four-pocket 
structure is essential for proper identification of relative 
magnetic ordering of spins of the two Fe atoms from the 
unit ceU in folded BZ Hi' 
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We assume that electron-electron interaction is short- 
range (Hubbard-like) and involves two couplings ~ be- 
tween fermionic densities from the same orbital and from 
different orbitals [37[. The Hamiltonian has the form 
H = H2 +^4, where 
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where where pi +P2 = Pa +P4, and Un is intra-orbital, and U12 inter-orbital interactions which we approximate 
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FIG. 2; Five relevant bare vertices. Solid and dashed lines 
represent fermions from c— band (near k = 0) and /—band 
(near k = (7r,7r)). 



by momcntuni- independent (on-site) values. 

The quadratic form can be easily diagonalized by 



'>Pi,p,a = cos9pCp^a + sin0pfp^a 
V'2,p,o- = cos6'p/p.cr - sin^pCp^o- 
with tan20p = 2rp/(e2,p — ei,p). This yields 

"2 ~ / ^ ^p'^p,a'^P,c + ^p/p,cr/p.o"; 
p,CT 



(2) 



(3) 



Under some conditions on the original dispersions ei,p 
and 62, p, and on the hybridization term Fp, the two 
bands of fermionic excitations form small hole and elec- 
tron pockets near (0,0) and Q — (tt, tt), with roughly 
equal size, as in the iron pnictidcs. This happens if, e.g., 
ei,p and e2,p change sign under p -^ p+Q, £1,0+62,0 > 0, 
and ./(ei.p-e2,p)2+4r2 < |ei_p + e2,p|. 
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Under these conditions, 



c ^ _ / 
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and ep describes a 



hole band with the maximum of energy at (0,0), while 
e^ describes an equivalent electron band with the mini- 
mum of energy at Q. Upon doping, chemical potential 
shifts, one Fermi surface gets larger while the other gets 
smaller, see Fig HI 

In itinerant systems, the interactions are expected to 
be small compared to the fermionic bandwidth, and the 
physics is dominated by fermions near the Fermi surface 
(FS). The projection of the Hubbard interaction term Hi 
onto c and / fermions leads to five different interactions: 
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where the momenta of c— fermions are near (0,0), the 
momenta of /—fermions are near (tt, tt), and the momen- 
tum conservation is assumed. We present these interac- 
tions graphically in Fig[2l 

We label the couplings with subindcx "0" to empha- 
size that these are the bare couplings. The terms with 



r(0) 



r(0) 



Ul"' and C/g are intraband interactions, the terms with 



r(0) 



r(o) 



Ui and U2 are interband interactions with momentum 
transfer and Q, respectively, and the term with U^ is 
interband pair hopping. 

Note that in our Fermi-liquid description, all vertices 
in Eq. ^ arc (5- functions in spin indices, i.e., all inter- 
actions are in the charge channel [38|, and there are no 
direct spin-spin interaction terms with spin matrices in 
the vertices. However, if the original Hubbard interaction 
is on-site, one can use another, equivalent, description in 
which Pauli principle is build into the Hamiltonian, and 
the intra-orbital terms with equal spin projections are 
eliminated from the Hamiltonian. In this description, 
C/f , C/4 , and C/g terms appear as effective Hubbard 



interactions, while C/j and C/3 appear as a magnetic, 
Hund term [H]. 

In explicit form, U^ are 



r(0) 



1 



Ur' = 2 [(t^ii + C^i2) - cos20oCOs20Q([/n - C/12)] 

U-^3 = ^^ (1 - cos 29q COS 26Iq) sin 26lo sin 26'q 

..(0) _ Uu + U12 Uii - U12 2 ^n 
U^ — ;:; \ Ti cos IUq, 



(0) _ t^ll + t^l2 Uii — c/12 



ur^ 



cos2 26'q. 



(6) 



For the case that we considered above (eip = 
-ei,p+Q,e2,p = -e2,p+Q,rp+Q = iFp), we have 9q = 
tt/2 T do, we have 



r(0) 



U^"> ^ U, 



(0) 



m 



(0) _ Uii + U12 Uii — C/12 2 



cos^26'o. 



C/f' = C/f ) = ^(1 + cos^ 20o) T ^ sin^ 29^ (7) 



We assume that the intra and inter-orbital Hubbard-type 
interactions Uu and C/12 are positive (repulsive). We see 
from d?]) that density-density couplings [/} ,U^ ' , and 



(0) 



5 are positive and the largest. The couplings C/j 



(0) 



and C/3 are smaller for the case when the hybridiza- 



tion term is even under p 



Q, i.e., Tp = r 



and are the same as Ul ,1/4 , and C/5 when Fp = 
Fp+Q. The first case corresponds to on-site hybridiza- 
tion and is more realistic that the second one, which 
requires hybridization to involve predominantly nearest 
neighbors. Below we will consider only the first case 
Fp-)_q. Note that in this situation, the sign of 



r(0) 



r(0) 



t/2 = C/3 depends on 6*0 and on the relative strength 
of the intra-orbital and inter-orbital Hubbard terms. If 
Uii > C/12 sin^ 26o / {1 + cos^29o) , these couplings are pos- 
itive, if C/11 < C/12 sin^ 26'o/(l + cos^29o), they are nega- 
tive. A more likely situation is when the intra-orbital 
Hubbard term C/n is larger than inter-orbital C/12, in 
which case C/j and C/3 are positive. 

For convenience, below we will be using dimensionless 
interactions 



U^No 



(8) 








FIG. 3: Diagrammatic representation of the equations for 
SDW, CDW, and SC instability temperatures. The equations 
'.'■•'J r,„A T^2w ^^^ obtained by adding and 

and for Ac 

and A*^„,, respectively. The equations for Tsc and Tsc ob- 
tained by adding and and subtracting equations for A*^ and 

A-''. 



for density-wave T^j^'J and T^ 



where Nq is the fermionic density of states (DOS) which 
we approximate by a constant. For itinerant systems, 
|u"| < 1 and can be treated within Fermi liquid theory. 
We will also count the momenta of /— fermions as devi- 
ations from Q (/p -^ /p+q) in which case all running 
momenta in the vertices arc small. 



introducing infinitesimal couplings 
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Density virave and pairing instabilities 

We searched for possible density-wave and Cooper- 
pairing instabilities for our model, and found that the 
ones which may potentially occur are spin density wave 
(SDW) and charge density wave (CDW) instabilities with 
momentum Q and with either real or imaginary order 
parameter, and superconducting (SC) instabihty either 
in pure s channel (the gaps Ac and A/ have the same 
sign), or in s"^ channel (the gaps Ac and Ay have op- 
posite sign). Density- wave instabilities with q = and 
pairing instabilities with q = Q do not occur within our 
model because the corresponding kernels vanish for a con- 
stant DOS. The instabilities with momentum-dependent 
order parameter, like a nematic instability [35[ also do 
not occur simply because we set all interactions to be 
momentum-independent and weak, and will neglect reg- 
ular (non-logarithmic) corrections which give rise to the 
momentum dependence of the scattering amplitude in a 
Fermi liquid [38|. 

The temperatures of potential density-wave and pair- 
ing instabilities are obtained by conventional means, by 



(9) 

with complex Agdw, ^cdw, and real A^/, and analyzing 
when the response functions diverge. We label the cor- 



n{r,i) rp{r,i) 



and 



responding instability temperatures as T^^'w '^cdw 

(s s^) 

Tsc ' , where r, i mean real or imaginary density- wave 
order parameter, and s, s+ mean s— wave or extended 
s— wave, respectively. The linearized equations for the or- 
der parameters are presented graphically in FiglS) They 
have non-zero solutions when 
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_rp{l\i) ST^ p(»-,'') 
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„(r.i) V^p(r,i) f 
cdw / J cdw I 



(10) 



Here 



sdw 



ui±U3, f|.^'^'' =1*1 =FW3 -2u2, 
U4 + U3, F^* ) = M4 - U3 



(11) 



are the full interactions in the SDW, CDW, and SC chan- 
nels. Eq. (flOl) is only valid for the largest instability 
temperature. Below such T, the ordering in one channel 
affects susceptibilities in the other channels. 
For the bare parameters as in ([7|) 
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M12 sin 26*0 

2ui2 sin^ 26*0 - mii(1 + cos^ 26*0), 
2 < 



2«!] 



P3 

— *— 



!] - 2m^ w ^12 sin^ 26*0 



+ M^«itii(l + cos^20o), 



W12 



sin^ 20n 
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(b) 



P3 



where un = UuNq, U12 ~ f/12-^0 are dimension- 
less intra-orbital and inter-orbital Hubbard couplings. 
The Stoner-like SDW and CDW instabilities require 
^sdw, ^cdw > 0. At the bare level, F^^^ is the largest 
positive interaction when wii(l + cos^ 26o) > M12 sin^ 26o 

(r) . 

and r^ijj^ is the largest positive interaction when uii(l -I- 
cos^ 29q) < U12 sin^ 29q i.e., the system undergoes a con- 
ventional SDW or CDW instability. The SC instabilities 
requires an attraction (a negative Tsc" ) and do not oc- 

( s) (s'^) 

cur at this level because both Tic and Tsc are positive. 



RG FLOW 

Beyond mean-field, the potential SDW and SC insta- 
bilities are determined by Ui at energies below the Fermi 
energy Ep, and generally differ from bare m° defined at 
energies comparable to the bandwidth, W. For small 
size of the FS, W >> Ep, and the intermediate range is 
quite large. At u^ < 1 the renormalization can be consid- 
ered in one-loop approximation. The one-loop diagrams, 
shown in FiglH contain particle-particle and particle-hole 
bubbles. The external momenta in these diagrams are 
of order running E > Ep, while internal momenta are 
generally of order W, i.e., much larger. In this situ- 
ation, the dependence on the directions of the external 
momenta is lost, i.e., a SC vertex with zero total momen- 
tum and an SDW vertex with transferred momentum Q 
are renormalized in the same way. The crucial element 
of our analysis is the observation that, for efj = —£„+„, 
particle-hole channel is undistinguishable from particle- 
particle channel, such that the renormalization in both 
channels are logarithmical and interfere with each other. 
The presence of the logarithms in both channels implies 
that the one-loop perturbation theory must be extended 
to one loop RG analysis for the running Ui (in the dia- 
grammatic language, one needs to sum up series of log- 
arithmically divergent parquet diagrams). The deriva- 
tion of the RG equations is straightforward (see Fig. [5]). 
Collecting combinatoric pre-factors for the diagrams, we 
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FIG. 4: Diagrams for one-loop vertex renormalizations. The 
renormalizations of ui and U3 are shown, others are obtained 
in a similar way. 



obtain 



Ui = U^ + M3 

M2 — 2u2{ui — U2) 

its = 2^3(2^1 - U2 - M4) 

M4 = — Uo — U4 



(13) 



where the derivatives are with respect to logW/E, and 
E is the running energy scale. Similar, though not iden- 
tical equations have been obtained in the weak-coupling 
studies of the cuprates with the "i— only" dispersion [39|. 
We see from Eq. p^ that the pair hopping term U3 
is not generated by other interactions, i.e., U3 = if 
U3 = 0. In the absence of U3, Tsdw and F^^" - Tgdw = ui 
increases and drives Tsdw up, while F*^'' = U4 logarith- 
mically decreases, as it is expected for a repulsive inter- 
action [40[. However, once M3 is finite, the system moves 
into the basin of attraction of another fixed point, at 




0.0 0.5 1.0 1.5 2.0 2.5 



4 


(a) ^ / 


saw 


2 


^^^■'' ,,'-' 


saw 




— «---„..--^"^^"^ 


p(s+) 





^^■^■^=----^- - 


sc 
sc 


-2 


\m 


cdw 


-4 


. \ 


^ cdw 



0.0 



0.2 



0.4 

L 



0.6 



0.8 



FIG. 5: (color online) The RG flow of Eqn. ^ in 
variables 114/113 and ui/u^. The fixed point is 114/113 = 
— 1/-\/5, iii/i'3 = l/\/5- The fourth variable, 112, becomes 
small near the fixed point compared to the other Ui, We used 
boundary conditions ul = 1*4, 1*3 = O.liti, and for simplicity 
set Wo = 0. 



which 



M3 oc 



1 



Hlnf 



«l = -li4 = ^,W2Cx|«3r/' (14) 



where u depends on the bare values of the couplings. 
In Figsin] and [H] we show the RG flow obtained by the 
numerical solution of Eq. ([1^)) . 

The two key features of the new fixed point are (i) \u^\ 
rapidly increases and eventually becomes larger than ui 
by a factor v5, and (ii) M4 decreases, passes through 
zero, changes sign, and then increases by magnitude and 
approaches —\uz\/V^ (see Fig. [5]) 

For positive Wg « Mj, these results imply that r^^.^, = 
Ml + 1*3 remains positive and the largest out of density- 
wave vertices i.e., the highest density- wave instability 
is a conventional SDW instability (see Fig. [5^). Note, 
however, that F^^^,^ is close second as it only differs by 
1*2 which under renormalization becomes relatively small 
compared to ui and U3 (u2 oc (1*3)^''^). The interaction 



1*4 — 1*3, becomes negative 



in the s+ SC channel, fI"/^ 

(attractive) below some scale (Fig. |BJi), while VW re- 
mains repulsive. We emphasize that the density-density 
vertex 114 changes sign under renormalization, becomes 
attractive and also supports SC. Moreover, the interac- 
tions in the SDW and the s+ SC channel F^^:^^^ = lii +1*3 
and Tsc = U3 — U4, become comparable to each other 
and eventually flow to the same value 1*3(1 -I- l/-\/5). The 
implication is that the SDW order and s+ superconduc- 
tivity are competing orders, determined by effective in- 
teractions of comparable strength. 

2\u2\ is the 
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strongest, positive, density-wave vertex, and F 
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liisi is a close second (see Fig. [SJd). F 



{s- 



is now positive (repulsive) , but F, 
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FIG. 6: (color online) The RG flow of the effective couplings 
in various density- wave and superconducting channels vs L = 
1*" log W/E. The boundary conditions are u? = 1*4 = 1*0, and 
1*2 = 1*3 = ±0.1i*o. The running couplings are in units of uq. 
Panel a) - the RG flow for 1*3 > 0. The extended s— wave 
(s+) superconducting channel becomes attractive above some 
L. The strongest density-wave instability is in SDW channel, 
for real order parameter. Panel b) - the same for 1*3 < 0. The 
conventional s-wave superconducting channel becomes attrac- 
tive above some L. The strongest density-wave instability is 
in CDW channel, again for real order parameter. 



sign under the renormalization and becomes negative (at- 
tractive), see Fig. ^. This implies that CDW now com- 
petes with a conventional s— wave SC. Near the fixed 
point, the interaction in the s— channel Vjc ~ ^l^sKl + 

I/V5) is now larger than T^^J^ ~ [""sKl ~ I/a/S) which 
implies that in this case s-wave SC likely wins over CDW. 
The generalization of this analysis to 4-band model 
(or even five) is straightforward and yields qualitatively 
similar behavior. 



Competing orders 



We next analyze in more detail Eqs. (fTO|) for 1*3 > 0. 
By construction, the upper limit of the integration over 
internal energies there is 0{Ep) as the contributions from 



higher energies are already absorbed into the renormal- 
ized vertices. When hole and electron Fermi surfaces are 
near-identical, i,e., e^ = — £^+0 holds down to the lowest 
energies, both SDW and SC susceptibilities are logarith- 
mic in T 

-TJ2 J de^Gi^^^Gi^^^^^ =TY,Jde^Gi^^G'Li,^_ 



tanh 



Lo \ dcu Ep 

— — = log 

2TJ uj ^ T 



intraband density-density interactions U4 and U5 become 
unequal). Second, solving the non- linear gap equation, 
we immediately find that the gap A obeys the same BCS 
relation 2A = S.SSTc ^^ for an isotropic s— wave state 
simply because the pairing kernel contains either two 
c— fermions or two d— fermions, but no cf pairs. 
„ The s+ and s SC states, however, differ qualitatively 
in the presence of non-magnetic impurities. For s— state, 
/-. suon- magnetic impurities do not affect Tc and non-linear 
gap equation 43| . For s'^ state, the impurity poten- 



and from pop the largest instability temperature is cither 



T, 



(r) 



sdw 



.(r) 



Epe 



rp(s 



Epe 



n( = +)| 



(16) 



(s+) 



As r^^(„ is still larger than T)tc ', the SDW instabihty 
comes first. This is what, we believe, happens at zero 
doping. Whether SC emerges as an extra order at a 
smaller T requires a separate analysis as the pairing sus- 
ceptibility changes in the presence of the SDW order. At 
a finite doping, all evidence is that the two FS become 
unequal, i.e., the condition e^ ~ ^^k+o breaks down. In 
this situation, the log 1/T behavior of the SDW polariza- 



(r) 

tion is cut, and T^^^ decreases and eventually becomes 

smaller than Tc . At larger dopings, Tc remains 
roughly doping independent, while magnetic correlations 
decrease. 

A remark about the SDW state. In the coordinate 
frame associated with folded BZ, Fe ions are located at 
I"! = inx,ny), where nx,ny are integers (we recall that 
we set interatomic spacing to one) , but also at r2 = {rix + 
1/2, Uy + 1/2). SDW instability with Q = (7r,7r) order 
antiferromagnetically spins within the sublattice where 
r = ri, and within the sublattice where where r = r2, 
but do not fix relative orientation between the spins in 
the two sublattices. To obtain full spin structure, we 
would need to analyze spin ordering within full four-band 
structure (two electron and and two hole orbitals), or go 
back into unfolded Brillouin zone. For localized spins, 
this type of order is described by Ji — J2 model for J2 > 
0.5 Ji. In the classical model, the angle between ri and r2 
sublattices is arbitrary, but quantum fluctuations select 
(0, tt) or (tt, 0) state [4l|,|4^. There is then an extra Ising 
degree of freedom, which was argued 35|, 1421 to remain 
broken even at T > Tsdw, when SU{2) spin symmetry is 
restored. 



SUPERCONDUCTING STATE 

The SC s+ state that we found has two features sim- 
ilar to a conventional isotropic s-wave state. First, the 
superconducting gaps on the hole and electron FS are op- 
posite in sign, but equal in magnitude. They, however, 
become unequal when Ep on the two FS become is differ- 
ent, which happens once the doping increases (or when 



tial [/i(q) has intra and intcrband components Ui{0) 
and Ui{'K), respectively. The Ui{'K) components scat- 
ter fermions with -|-A and —A and acts as a "magnetic 
impurity" [4J, |45| . Specifically, for the s"*" state, normal 
and anomalous Greens functions in the presence of im- 
purities are 



G\ 



■cj 



Zuj^UJm ± Ck 



F, 






±- 






'ziAco?, 



J + 4 



(17) 



and the fermionic Z ^ I + T,{ujm)/ujm and the renor- 
malized gap A^ in the Born approximation are given 
by 



Z ^ 1 

A, 



C/,:(0) + [/,(7r) 



VA^ 
-1 = - 



for A, 



A2 



(18) 



where A = A(T) is the frequency-independent order pa- 
rameter, and &T = 2t/i(7r)/A(T). Below we use &t=o = b 
as a measure of the strength of impurity scattering. Note 
that 6 is a complex function of the impurity strength as 
the order parameter is also affected by impurities (see 
below) . 

For C/i(7r) = 0, A = A, i.e., superconductivity is not 
influenced by impurities. For Uiin) 7^ 0, Ai^^ becomes 
frequency dependent, as if the impurities were magnetic. 
At r = 0, and b > 1, the system displays gapless 
superconductivity [46[: in real frequencies A^^^ ex —iuj 
at small to, and the DOS at zero energy acquires a finite 
value N{lu = 0) = (1 - (1/6)^)^/^ Superconductivity at 
T = eventually disappears when A vanishes, i.e., when 
b tends to infinity. 

The parameter A can be re-expressed in terms of 
Aq(T), which is the BCS gap in the absence of impuri- 
ties, and 60 = 2C/i(7r)/Ao(T), which hnearly depends on 
the impurity strength. The relation between S = A/Ao 
and bo (and between b = bo/5 and 60) is obtained from 
the self-consistent condition on the order parameter 



^e// 



A2 
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FIG. 7: Tc (a) and bo = 2U^{^v)/Ao (b) as functions of 
b = 2[/i(7r)/A, where A is the order parameter, and Ao is the 
gap in the absence of impurities. The inset in (a) shows the 
dependence of Tc on bo . 



An = 



'e// 



Ao 



VAT 



(19) 



where u^^ff is the normaUzed interaction in the s"*" chan- 
nel. Solving these equations, we obtain after some alge- 
bra the relations which express bo in terms of b. They 
are 



6o = bexp [ — — 



b< 1 



ba ~b cxp 




— sin ^ cosh fe , 6 > 1 

2 b J 

(20) 

The first regime corresponds to S = b^/b > 5o, and 
holds for 6o < e~'^/^ « 0.465. The second regime cor- 
responds to 6 > 1 and describes a gapless superconduc- 
tivity (A^^^o = 0). When the order parameter A tends 
to zero, and b tends to infinity, b^ approaches 0.5. We 
plot 6o vs b in Fig. [Tfb). 

The vanishing of superconductivity at T = when 6o 
approaches 1/2 also follows from the generic dependence 
of Tc on the impurity strength. The calculation parallels 
the one for an s— wave superconductor with magnetic im- 
purities [46| and yields 



Tr V 2 All Tr 



-i^[\ 



(21) 



where ip{x) is the diGamma function. One can easily 
check that Tc vanishes when 6o approaches 1/2. We plot 
Tc{b) and r,(&o) in Fig. [Tta). 



Spin response of an s^ superconductor 

The dynamical spin susceptibility of a superconductor 
is given by an RPA-type formula 



Xs(q,^) 



x2(q,^) 



^{r) 



l-mx°(q,^) 



(22) 



where Xs('Zi ^) is the (dimensionlcss) susceptibility of an 
ideal s"*" SC (the sum of GG and FF terms with spin 



matrices in the vertices) 



In our case, when e^ 






and the gap changes sign between hole and electron FS, 
one can easily verify that Xs(q ~ Qj^) coincides with 
the particle-particle susceptibility for either c— or /— 
fermions. This leads to several consequences. 

1. In the normal state, Xs(Qi^) = log^F/(— «f^), 
that is ImXs(Q,tt') only weakly (logarithmically) 
depends on frequency. This could be verified in 
INS experiments. 

2. In a superconducting state, Xs(Qj^) has a reso- 
nance below 2A. Indeed, at T = 0, in the clean 
limit and at small 57 and q — Q, 



X°(q,r!)=log 



El 

En 



1 



4A2 



{n^-v^iq~Qf){23) 



where v = vp/V^ is the velocity of the Anderson- 
Bogolyubov mode in two dimensions (2D), and Eq 
is the largest of A and the cutoff energy associated 
with non-equivalence of the two FS. Substituting 
this into Xsiq, ^), and assuming F^^^ \ogEp/Eo < 
1, i.e., no SDW instability, we find the resonance at 

n - V("'(q-Q)' + f^o. where flo = 2A(l/ril- 
logEp/EoY^^. This resonance has been earlier ob- 
tained in the numerical analysis in Refs. [23|, |47|. 
It bears both similarities and differences with the 
spin resonance in dx2_y2 SC. On one hand, both 
are excitonic resonances, and both occur because 
the gap changes sign between the FS points k and 
k -I- Q. On the other hand, the resonance frequency 
in a d^2_y2 SC disperses downwards because of the 
nodes, while for a nodeless s^ SC, the resonance 
disperses upwards, with large velocity. Indeed, this 
is only valid if fJ ^ 2A, otherwise the dispersion 
becomes more complex. 

Note in passing that, because the two gaps have 
opposite signs, there should also exist a resonance 
mode in the particle-particle channel, at momen- 
tum k = Q, similar to the Leggett mode in a two- 
band superconductor [48|. 

3. An s"*" superconductor has a rather peculiar low- 



frequency behavior of Imxs (q ' 
clean limit, 



Q,fi^O). In the 



Imx°(q,f7) 



n 



ex 

0=0 k 



E^^^ 



' dE^ 



(24) 



k. ME) 

'^k'^k+q 



is Fermi func- 

A^Af)/{E^Ei^^) 



where E^ ~ V^A^ + ef. 
tion, and Ck,q ~ ^ + {^i 

is the coherence factor. We see that the coher- 
ence factor vanishes identically for q = Q such that 
Imx°(q,fi)/O|o=ocx(q-Q)2. 



NMR spin-lattice relaxation rate 



The spin-lattice relaxation rate measured by nuclear 
magnetic resonance (NMR) is given by 



^-^E 



ImXs(q,^) 



Ti 



9. 



rE^'(q'O) 



0=0 




Imx^(q,f]) 



9 



(25) 



0=0 



Because XsCQj^ = 0) is enhanced near Q, this region 
contributes most to the momentum sum. The smallness 
of ImXs(q, ^)/^|n=o for q~ Q then implies that l/Ti 
has extra smallness in a clean s+ SC [by the same reason, 
there is no Hebel-Slichter peak in l/Ti near T^. 

In the presence of impurities, 

ImXs(Q, f^)/i^|n=o remains nonzero, and l/Ti « 
rimXs(Q,^)/fl|f2=o/c^^'?Xs(q,^ = 0). The fuU expres- 
sion for l/Ti is rather involved as one has to include 
the full G and F, and the full vertex. It simplifies 
considerably if we neglect vertex corrections and assume 
that intraband scattering C/i(0) (harmless for supercon- 
ductivity) well exceeds A. In this case, wc obtained 
analytically, at a finite T, 



1 



1 



dx 



lAI 



4 cosh 



2T 



V(|A|2-x2)2+4x2(A")2 



(26) 



where A is given by Eq. (fT5|) with BCS T— dependent 
A(T). We verified numerically that lowest order vertex 
corrections do not change the result in any significant 
way. 

In FigO we plot the normalized temperature depen- 
dence of l/Ti(T) for several values of 6 = 2f/,(7r)/A(r = 
0). Stronger impurity scattering corresponds to larger 
b [it doesn't make a difference wthether to parametrize 
the impurity strength in terms of h, which depends on 
impurity strength in a complex way, or &01 which scales 
linearly with the impurity strength, because of one-to- 
one correspondence between h and 60, see Eq. (|20p ]. 

For 6 < 1, the low-T behavior is exponential, as is ex- 
pected for a superconductor without nodes. However, we 
see that for h > 0.3, there is a wide intermediate T range 
where the behavior of l/Ti closely resembles a power- 
law T". The exponent a decreases as b increases from 
a « 3 for 6 = 0.3 to a « 2 for 6 = 0.9. The T^ behavior 
was suggested based on experimental fits and was pre- 
sented as evidence for d— wave superconductivity in Fe- 
pnictides. Our results show that l/Ti(T) in a dirty s+ 
superconductor mimics a power-law over a wide T range 
even when the DOS still vanishes at w = 0, and Tc is only 
slightly affected by impurities. Furthermore, we argue. 




0.4 0.6 0.8 1 
T/T 




FIG. 8: (color online) (a) Calculated temperature dependence 
of l/Ti for an s^ superconductor with nonmagnetic impu- 
rities. The normalization is chosen such that l/Ti = 1 at 
T = Tc- The theoretical curves are for various values of the 
parameter b = 2Ui(ji)/ /S.{T = 0) which measures the strength 
of the pair-breaking component of non-magnetic impurities. 
Gapless superconductivity occurs for 6 > 1. The experimen- 
tal data are taken from Refs. [33. |34J . (b) Theoretical 1/Ti 
for different b vs power-law forms . All theoretical dependen- 
cies are exponential in T at very low T, but are described by 
power-laws T" over a wide T range below Tc- The exponent 
a decreases as 6 increases from a ~ 3 for 6 = 0.3toa~2 for 
6 = 0.9. 



based on Fig. [HKa) that the experimental T dependence 
of I/Ti can only approximately be fitted by a particular 
power of T. We believe that the reported power-law form 
reflects intermediate asymptotics of a complex T behav- 
ior of l/Ti, and one should reduce temperature further 
to be able to distinguish between a true power-law and 
exponential behavior [49|. 

Note in passing that the theoretical behavior is expo- 
nential at the lowest T only if 6 < f . For & = 1, which 
is the critical b for a gapless s+ SC, 1/Ti ex T^^/^ at the 
lowest T, and for larger &, l/Ti(T) ex T. 
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1.0F 




Uniform susceptibility 



FIG. 9: (color online) Calculated temperature dependence of 
the uniform susceptibility for various values of b. The exper- 
imental data are taken from Ref . [3411 . 



Finally, we also computed uniform spin susceptibility 
Xs{T) « Xs(9 = ^ = 0)i measured by Knight shift. It is 
obtained by standard means [50(] , and for a superconduc- 
tor with s^ gap symmetry is given by 



XsiT) 



Xs{Tc) 
Xs{Tc) 



1- 
1- 



tanh -— - Im- 



2T (A2-2;2)3/2_2t/,(7r)x2 



tanh ■ 



Ini- 



2T A(T) (1 - ^2)3/2 „ 5^ 



(27) 



where Ux = 2;/A and A depends on x and is given bv (|18p 
We emphasize that ladder series of vertex corrections 
must be included in the calculation of Xs{T) to recover 
SU{2) spin symmetry. Observe that Ui{0) drops from 
the expression for Xs(T), because impurity scattering of 
electron within cither hole or electron FS does not dif- 
ferentiate between a conventional s—wavc and s'^ gap 
symmetry. 

Eq. ([27]) is similar, but not identical to the expres- 
sion for Xs (T) in an ordinary s-wave superconductor with 
magnetic impurities [50|. In both cases, Xs{T) differs 
from free-fermion value. However, for magnetic impuri- 
ties, Xs(T = 0) becomes finite for any nonzero strength 
of the impurity scattering, while in our case, the impuri- 
ties are actually non-magnetic, and XsiT) still vanishes 
at T = for all 6 < 1, for which the DOS still vanishes 
at zero frequency. 

We plot Eq. (gT]) in Fig. [5] for the same b as 1/Ti. 
We see the same trend as in Fig. [8j the theoretical T 
dependence of Xs{T) is exponential in A/T at the lowest 
T, but rapidly deviates from exponent already at small 
T, and is roughly a power-law in T in the same T range 
where l/ri(T) can be fitted by a power-law. This is 
another indication that one should perform Knight shift 
and 1/Ti measurements down very low T to be able to 
distinguish between the nodeless s"*" state and a SC state 
with gap nodes. Note also that the same 6 = 0.7 which 
fits 1/Ti data by Matano et al [3J| also fits reasonably 
well their data on the Knight shift. 



CONCLUSIONS 

To conclude, in this paper we presented Fermi liquid 
analysis of SDW magnetism and superconductivity in 
Fe— pnictides. We considered a two-band model with 
small hole and electron pockets located near (0, 0) and 
Q = (tt, tt) in the folded BZ. We argued that for such 
geometry, particle-hole and particle-particle channels are 
nearly identical, and the interactions logarithmically in- 
crease at low energies. We found that the interactions in 
the SDW and extended s— wave channels (Ak = — Ak+q) 
become comparable in strength due to the increase of 
the intraband pair hopping term and the reduction of 
the Hubbard-type intraband repulsive interaction. We 
argued that at zero doping, SDW instability comes first, 
but at a finite doping, s"*" superconducting instability oc- 
curs at a higher T. 

This s"*" pairing bears similarity to magnetically me- 
diated dj;2-y2 pairing in systems with large FS with hot 
spots in the sense that in both cases the pairing comes 
from repulsive interaction, peaked at Q, and requires the 
gap to change its sign under k ^ k + Q. the difference 
is that for small pockets, the gap changes sign away from 
the FS and remains constant along the FS. 

We analyzed spin response of a clean and dirty s"^ 
superconductor and found that (i) it possess a reso- 
nance mode which disperses with the same velocity as 
Anderson-Bogolyubov mode, (ii) intraband scattering by 
non-magnetic impurities is harmless, but intcrband scat- 
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tcring affects the system in the same way as magnetic 
impurities in an s— wave SC, (iii) 1/Ti has an extra small- 
ness in the clean hmit due to vanishing of the coherence 
factor, (iv) in the presence of impurities, there exists a 
wide range of T where the T— dependencies of 1/Ti and 
the uniform susceptibihty for an s+ SC resemble the ones 
for a SC with nodes. 

Note added While completing this work we became 
aware that similar results for spin-lattice relaxation rate, 
1/Ti, in the superconducting state have been obtained in 
Ref. 0]. 

We acknowledge helpful discussions with Ar. Abanov, 
J. Betouras, O. Dolgov, K. Ishida, M. Korshunov, I 
Mazin, J. Schmalian, Z. Tesanovic, and G.-q. Zheng. 
Ave has been supported by NSF-DMR 0604406 and MPI 
for Physics of Complex Systems in Dresden. 
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